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Abstract 

^ . Greedy Maximal Scheduling (GMS) is an attractive low-complexity scheme for scheduling in wireless networks. 

Recent work has characterized its throughput for the case when there is no fading/channel variations. This paper 



aims to understand the effect of channel variations on the relative throughput performance of GMS vis-a-vis that 
of an optimal scheduler facing the same fading. The effect is not a-priori obvious because, on the one hand, fading 
could help by decoupling/precluding global states that lead to poor GMS performance, while on the other hand 
fading adds another degree of freedom in which an event unfavorable to GMS could occur. 

We show that both these situations can occur when fading is adversarial. In particular, we first define the notion 



en 

of a Fading Local Pooling factor (F-LPF), and show that it exactly characterizes the throughput of GMS in this 



setting. We also derive general upper and lower bounds on F-LPF. Using these bounds, we provide two example 
networks - one where the relative performance of GMS is worse than if there were no fading, and one where it is 



H ■ better. 

Index Terms 

Local Pooling factor, Greedy Maximal Scheduling, Throughput Region, Channel Fading. 

I. Introduction 

This paper analytically investigates the effect of fading on the throughput performance of a natural and popular 
scheduling algorithm: Greedy Maximal Scheduling (GMS) |[T2l . El, (H, |5j. As with any scheduling algorithm, 
GMS is a way to determine which wireless links can transmit at any given time, based on their mutual interference 
characteristics and their current level of fading. In particular, GMS involves first associating a weight with each 
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Proceedings of IEEE Infocom, Orlando, FL, March 2012 1131 . 



link - which depends on the load of the link and its channel condition. Then, GMS involves iteratively turning on 
the heaviest link that does not interfere with links already turned on. This is repeated every time slot. 

GMS has empirically shown to have very good throughput and delay performance; recent theoretical advances 
HI, 0, lfT5l , HI, Q, ifTOl characterize its throughput. All of these works assume that there is no fading; ie that 
the rate a link can support is invariant as long as all the links that interfere with it are not simultaneously on. Our 
work investigates what happens to this performance in the more realistic setting with intrinsic channel fading as 
well. In particular, we compare the relative throughput of GMS as compared to that of an optimal scheduler. 

Our results demonstrate that the effect of fading is quite subtle; in particular, in some instances fading can 
degrade the relative performance of GMS, while in other cases it can improve it. The former reflects the fact that 
fading provides an extra degree of freedom and complexity in the system, which GMS may not be able to handle 
as well as in a system without this fading. The latter reflects the, perhaps more subtle, fact that the sub-optimality 
of GMS (even without fading) is tied to the existence of special global system configurations that result in poor 
performance. The presence of fading "breaks up" these global configurations - not allowing them to occur too 
often - allowing GMS to perform relatively better. 

Specifically, our contributions are as follows: For a given wireless network with fading channels, 

1) We define a new quantity, called Fading-Local Pooling Factor (F-LPF), analogous to LPF defined in (H 
that characterizes the performance of Greedy Maximal Scheduling (GMS) in wireless networks with fading 
channels. Furthermore, we show that Fading-LPF is a lower bound on the fraction of throughput that can be 
stabilizable by the GMS when the arrivals and channels are independent and identically distributed over time. 

2) With arbitrary arrival and channel state process, we show that Fading-LPF is an upper bound on the fraction 
of throughput that can be stabilizable by the greedy schedule. More specifically, we construct an adversarial 
arrival and channel process with long term averages that lie outside the scaled throughput region and show 
that GMS policy cannot stabilize the queues. 

3) We further provide lower and upper bounds on Fading-LPF that are easy to evaluate. We provide two example 
networks with specific fading structure and use the derived bounds to demonstrate that fading can either 
enhance or degrade the relative performance of GMS as compared to the non-fading scenario. 

4) With fading, we can represent the channel model as a collection of global channel-states, where each state 
is associated with an independent set and an occurance probability. A natural question that arises is the 
following: Is the acheivable rate-region with fading simply the (channel-probability weighted) average of the 
per-state scaled rate regions, with the scaling parameter simply being the conventional LPF for each state? 
We show that this is in general not true. However, we derive a region that can be stabilized by the GMS in 
wireless networks with fading channels. This region is characterized based on the interference degree of the 
subgraphs (generated from original network) and the fading distribution. 



A. Related Work: 

Transmission scheduling has been a key challenge in modern wireless systems. The MaxWeight algorithm, 
proposed in |[T4l . has been the inspiration for many approaches to address this in various wireless systems (see |@] 
for several variants). However, this algorithm suffers from centralization as well as computational complexity. 

Thus, there has been significant research in finding sub-optimal (i.e., achieving a subset of the throughput region) 
distributed scheduling algorithms with low complexity. The authors in |fl2 1 propose one such policy called Greedy 
Maximal Scheduling, whose time complexity is linear in the number of links, and has a distributed implementation 
(8). There are other sub-optimal, randomized algorithms that have been proposed with similar performance as GMS 

ED, 0. 

The authors in Q have been the first to study the performance of GMS under a general interference model. 
They have identified conditions (so called 'Local Pooling') under which there is no loss in the network throughput 
region with GMS. The notion of Local Pooling has been extended to a multi-hop regime by 03J. 

This condition being identified as too restrictive, the authors in have defined a new quantity called Local 
Pooling Factor (LPF) that exactly characterizes the fraction of throughput region achieved by GMS, and show that 
over tree networks with a K— hop model for interference, GMS achieves the entire throughput region. Additional 
characterizations, including a per-link LPF |]9] and bounds to characterize the stability region ITTOl . have been 
proposed in literature. 

The authors in [T| exactly characterize, using graph theoretic methods, the set of network graphs (with only the 
primary interference constraints) where GMS is optimal (LPF = 1). Finally, the authors in Q have studied the 
performance of GMS with the SINR interference model, and have shown that GMS exhibits zero LPF in the worst 
case. 

All the above results assume that there are no channel variations (fading). In this paper, we study the effect of 
channel variation on the performance of GMS. 

II. System Model and Back Ground 

We consider a wireless network consisting of K links labeled as {1, 2, 3, ..., K}. Let K denote the set of links 
in the network. Each link I consists of a transmitter and receiver. We assume time to be slotted. Each time slot is 
composed of two parts. The first (control) part is reserved for making the transmission decision and second part 
for transmitting the packet. At time slot t, we denote the channel capacity of link by Cr\t]. We assume that the 
capacity varies from slot to slot, and is constant during a time slot. We consider collision interference/protocol 
model and denote the set of links that interfere with link I by X[. We say that the transmission on link I at time t 
is successful, if no link in the xi transmits during the same time t. The maximum number of packets that can be 
successfully transmitted in time slot t on link I is bounded by Ci[t]. 

We assume single hop flows in the network. Let A\ [t] denote the number of packets that arrive at transmitter of 
link I at time slot t. We assume that arrival processes is bounded and average rate of arrivals for link I is denoted 



by A/. 

For simplicity we first consider ON/OFF channels (i.e Ci[t] = Oorl) and later show that our results can be 
extended to channels with finite number of channel states. For the ON/OFF setting, global state (GS) refers to 
specifying the set of links that are in 'ON' state. Let GS(t) denote the set of links that are in 'ON' state in time 
slot t. Let it (J) denote the fraction of time the network is in global channel state J, where links in set J are 'ON' 
and links in the set K.\J are in 'OFF' state. Let tv :- {it (J), J C /C} denote the fading structure. 

Assumptions. : 

Al (Long-term Averages): We assume that the long-term time averages of arrivals and channel states satisfy the 

following: 

1 T 

-J2Mt\^^i as T^oo. (1) 

and 

1 T 

J, Y^ 1 GS(t)=J -> n(J) as T ->• oo. (2) 

i=0 

A2 (Randomness): We assume that arrivals are mutually independent Ltd processes with Xi = E[A[[t]]. Similarly 
the channels are independent across time and form a stationary process with it (J) = E[]-GS(t)=j]- 

While both assumptions Al and A2 specify the same long-term averages, we note that assumptions in Al allow 
for arrival and channel state processes to be dependent across time and across links in a deterministic, and possibly 
adversarial manner. The necessity for the above sets of assumptions will be clear as we state our main results in 
Section HH 

A. Preliminaries 

As discussed earlier, there is a rich history of analysis of GMS algorithms for the non-fading case (3), 0, J9], 
ifTOl . iPTI . 0. In this section we build on this notation in literature to allow for time-varying (fading) channels. 

We define Interference graph XQ for a set of links as follows: Each link is represented by a node and an edge 
is drawn between two nodes if transmissions on the corresponding links in the original graph interfere with each 
other. This model captures many existing wireless models and is quite general. We define the Independent set on 
this graph as set of nodes with no edges between them. Let Q\ [t] denote the number of packets present at the 
transmitter at time t waiting to get scheduled on link /. Let Si[t] € {0, 1} denote the schedule decision for link I at 
time t. At each time t, a schedule S[t] is determined based on the global queue state and channel state information 
at time t, that is (Q[t]), C[i\). We also assume that arrivals occur at the end of time slot, thus we have the following 
queue dynamics: 

Qi[t + 1] = (Qi[t] - C / [i]5,[t]) + + Mt], (3) 

where a + = max(0, a). 



Given the arrival traffic rate {Xi}i^c ar >d a scheduling policy, we say that the network is stable under scheduling 
policy if the mean of the sum of queue lengths is bounded. We say that an arrival rate vector {Xi}i^c is supportable 
if there exists any scheduling policy that can make the network stable. We call the set of all arrival vectors that 
are supportable by throughput region and denote it as Aj, where / denotes that the channels are fading. 

We say that a scheduling policy is throughput optimal if it can stabilize the network for all arrival rates inside 
the throughput region. 

Definition 1: ([5]) The interference degree dr(l) of link / is the maximum number of links in the set {I UZ/}fhat 
can be active at the same time with out interfering with each other. The interference degree di(G) of a graph 
G = {V, E} is the maximum interference degree across all its links in E 

Consider a wireless system with 4 links. Let Z\ = {2}, T-i = {1,3,4}, Z3 = {2,4} and Z4 = {2,3}. The 
interference graph is shown in the Figure Q] with the corresponding dj(l). The interference degree of this example 
graph is 2. 



O 1 




d/1) = 1 

dj(2) = 2 

d^) = 1 



Fig. 1. Interference Graph where nodes denote the links and edges denote the interference constraints. 



Definition 2: Given an interference graph, an independent set corresponds to set of nodes (links in the original 
graph) such that there is no edge between any two nodes in the set (no two links interfere in the original graph). 
Further, it is maximal if it is not a subset of any other independent set. For a set of links L, define a matrix Ml 
whose columns represent the maximal independent sets on the set L, with \L\ rows one for each link. We assume 
links are naturally ordered and rows in Ml are assigned according to the defined order. For J C L, let Mj^l denote 
the matrix with \L\ rows and is constructed from Mj as follows: columns from Mj are used and zero row vectors 
are added for links which do not belong to set J. Let CH(Mj^l) denote the convex hull of all column vectors of 
matrix Mj^l- 

For the above example with 4 links, let J = {1, 2, 3} and L = {1, 2, 3, 4}, we have 

(1 o\ 



Mj 



V 



1 

1 



and 

(l o\ 



Mj, 



1 

1 

Note that the set A^ := {A : A < /2; fl G CT-L(Ml)} characterizes the throughput region of set of L links if no 
fading were present. We now define the throughput region with the fading structure, 

Definition 3: The throughput region Aj for a given network with fading pattern vr( J) is described as follows, 

A/ = {X : A > 0, A < J2 vr(J)r/j where fjj G Cft(M JjX :)}. 

J 
Definition 4: ((51) The efficiency ratio j* ol under a given scheduling policy is defined as follows, 

7* oi = sup<^ 7 : the policy can stabilize for allthe arrival rate vectors A G 7 Ay >. 

Definition 5: Given x(J) G [0, 1], we define a new region Af(x) as follows, 

A / (f) = |X : A > 0, A < Y^ -x ( J) vr ( J) r/j where ffj G CH(Mj jX; )|. 

J 
Note that throughput region is same as Aj(l). 

B. GMS Algorithm SJTjl 

We now describe the Greedy Maximal Scheduling(GMS) Algorithm. GMS essentially finds a maximal schedule 
in a greedy fashion. Each node in the interference graph is assigned weight equal to f(Qi(t)Ci(t)), where /(.) 
is a strictly increasing function that is zero at and tends to infinity as Qi(t)Ci(t) — > 00. It then proceeds as 
follows: it finds the node with maximum weight in the whole network and adds it to GMS schedule (ties are 
broken arbitrarily), it further discards all the neighboring nodes along with the selected node and repeats the above 
procedure on the reduced graph, till there are no more nodes left in the interference graph. 

III. Main Results 

In this paper, we characterize the performance of GMS algorithm for wireless networks with time-varying 
channels. We define the fading local pooling factor, a* L (iv), for a set of links L(C /C) with fading structure n 
as follows: 

g* l {-k) = inf{o- : 3 (f 1 , <f 2 G $(L) such that acfi > <f 2 }, (4) 

where, 

$(L) = {$:$= J2 < J )VJ where ffj G CH(M JnLiL )}, (5) 

J-.JCK. 
and Fading-Local Pooling Factor (F-LPF) for a network G, 0-^(71"), with fading structure iv as follows: 

ctg(tt) =min L:L c/ccr2(7r), (6) 



Note that the above definition reduces to the known definition of LPF for a graph ||5) when there is no fading, 
i.e, when vr(/C) = 1. 

The F-LPF can be understood as follows: Consider arrivals only to links of set L (assume arrivals to other links 
are 0); when the links in set J are 'ON' (others are 'OFF'), GMS will pick a maximal schedule among the 'ON' 
links, i.e. a column of Mjc\L,l- Thus vector ffj is the long run average of these maximal schedules when system 
is in state J; so fjj € C%{Mj^l,l)- Thus &(L) is the set of all long-run average service vectors that could appear 
due to GMS when the arrivals are restricted only to set of links in L. For any two vectors <pi,(j)2 £ &(L>), it may 
thus happen that GMS results in </>2 service vector, when it should have been (j)\ (for the optimal case). Thus o~ l {tt) 
is the worst possible ratio difference among all the possible service vectors of $(X). 

Dual Characterization and Implications: In the same spirit as (3], 0, the Fading- Local Pooling Factor has a 
dual characterization, as noted in Lemma |3l and displayed below. The F-LPF, cf l (it), is given by the solution to 
the following optimization problem: 

a* L (ir) = max V vr L (J)a(J) (7) 

x,a(J)MJ) j^ L 

s.t : x'M, lL > a(J)e VJCL 

x'M JiL < b{J)e' VJCL 

Y^ vr L (J)6(J) = 1, (8) 

J-.JCL 

where e is a column vector of all ones, (•)' is the vector transposition operation and 7Vl denotes the marginal 
distribution on set of links L induced by 7r. 

Observe that each fading state J induces a network defined by ON edges (i.e., all OFF links are removed from 
the network). Thus, one could ask if with fading channels, the F-LPF can be determined simply by computing 
the "standard" LPF (denoted by c*(J)) for each of these induced networks, and then averaging these quantities 
(weighted by the steady-state fractions of times for each of the fading states) over all possible fading states? In 
other words, is the following true? 

*!(*) = £ n L (J)a*(J) 

J-.JCL 

where a* (J) is the standard LPF Q for the network that is induced by state J. 

An important insight that emerges from the dual characterization is that such averaging does necessarily not 
hold, in particular because the possibly adversarial nature of the fading channel does not permit averaging. Note 
that the adversary cannot change the long-term fractions of the global states - it can merely change the temporal 
correlations. Inspite of this, averaging does not hold, as clearly shown in Example B in Section IIII-BI) . 

In a tree network with fading as in Example B (see Section ITlI-BI ). while the LPF for each state is '1', the F-LPF 
is less than 4/5 which is lower than any convex averaging of the states! This discussion implies that the regular 



LPF does not immediately extend to the case with fading. This motivates us to explicitly develop the local pooling 

factor in the presence fading, and understand its implications. 

Contributions: 

A. Characterization in terms of F-LPF: 

Our first contribution, Theorem [Q characterizes the efficiency ratio of GMS algorithm in the presence of fading. 

Theorem 1. a) (Upper Bound) Under a given network topology and channel state distribution with Assumption 
Al on the arrivals and fading channels, the efficiency ratio of GMS (j*) is less than or equal to (Jq{tt). 

h) (Achiev ability) Under a given network topology and channel state distribution tv with Assumption A2 on the 
arrivals and fading channels, the efficiency ratio of GMS (7*) is greater than or equal to o~q(iv). 

Implications: The above result enables us to understand the performance of GMS compared to the optimal 
scheduler in the presence of fading. In particular, computing bounds on <Tq(it) leads to insights on the positive 
and negative aspects of fading (discussed further in Theorems [2] and [3]). Observe first that as long as the long-term 
averages on the arrivals and channels are satisfied (Assumption Al), we can construct an arrival and channel process 
that ensures that the efficiency cannot exceed the F-LPF &q(it). Further, for typical arrival and channel processes 
with sufficient randomness (in this paper i.i.d. assumptions have been imposed, however this can be weakened), 
the converse holds wherein o-q(tv) is achievable. 

Proof Discussion: For the first part, we extend the ideas in [0, to construct an adversarial arrival and fading 
process pattern when arrival rates are outside the (ctq(7t) + e)Af and show that a set of queues are unstable 
under GMS policy. For the second part, we use the approach in (3), as follows: we show that if A is inside 
(q-q(tv) — e)Af then GMS policy can stabilize all the queues in the network. We look at the deterministic fluid 
limit of the system and exhibit a Lyapunov function whose drift is negative under the GMS policy. We have that 
fluid model is stable and therefore that the original system is stable. 

Theorem 2 (Upper Bound). For every J C K and any (jlj, vj, Hj) such that jlj, vj E CH(Mj), vj < Hjjlj, we 

have that 

., w Zjac*(J)HjHj(l) 
a G 7T < maxi — ^ — — -, 

where fXj(l) = if I ^ J. 

Implications: While (Tq(tt) is defined only though an optimization problem, the upper bound permits an explicit 
solution. This bound is useful, as evidenced in Example B provided in Section ITlI-Bi In particular this upper bound 
is useful to illustrate that the F-LPF is not a simple convex combination of the standard LPF averaged over the 
fading states, and that adversarial fading can indeed worsen the performance of GMS. 

Proof Discussion: Though the proof follows from straightforward algebraic computations, the value of the theorem 
lies in the smart selection of (jlj, vj, Hj) vectors that satisfy the inequality stated in the above theorem. In the 



worst case the bound yields 1; however we can use the existing results in literature |H] to get good bounds. Thus, 
the tightness of the upper bound depend up on the ability to identify good vectors that satisfy the above constraints. 

Theorem 3 (Lower Bound). 

aL ^' ~ V " nr (T\AT(M V (9) 

l^JCL 7T L{ J )N(Mj) 

where n(M) = mm j ^2- Mi j , N (M) = max,- Ylj Mjj. ttl denotes the marginal distribution on set of links L 
induced by iv and can be computed as follows, 



tt l (j) = Yl »co 



I:ICK,InL=J 

Implications: The ability to compute a lower bound leads to the interesting observation that fading can help 
improve efficiency. This is because, by turning links 'OFF', fading "breaks up" some of the bad global states that 
can lead to poor GMS performance. This is explicitly brought out in Example A in the context of a six-link network. 

Proof Discussion: The lower bound is derived using the dual formulation of the F-LPF, see (|7]). We find a point 
in the dual search space that satisfies all the constraints in the dual characterization, thus yielding a lower bound 
on the primal problem. Observe that n(Mj) corresponds to the minimum number of links that needs to be 'ON' 
in any maximal schedule on set of J links and N(Mj) denotes the maximum number of links that could be 'ON' 
among all the maximal schedules on set of J links. Thus, the lower bound can be computed easily and can be 
shown to be tight for some wireless networks. As an interesting aside, note that the lower bound provided is always 
better than the inverse of the interference degree of graph G (see Corollary 1). 

We now present two examples: A and B, one in which fading reduces the relative performance of GMS and the 
other in which fading enhances the relative performance of GMS respectively to illustrate the value of the above 
results. 

B. Examples: Benefit and Detriment with Fading 

Example A: A network where fading structure improves the relative performance of GMS: Consider a graph 
with six links K = {a, b, c, d, e, /}. The interference graph for the six links is shown in the Figure IIII-Bi Each 
link is either is state 'ON' or 'OFF'. We consider the following fading structure, 7v, for J C K. 

7T(J)=pl J l(l- P ) 6 -| J l, 

where \J\ denotes the size of set J. Note that p = 1 corresponds to the no-fading case. 

Using our results, we compute the lower bound and upper bounds on local pooling factor cr£.(7r) and is plotted 
in Figure IIII-Bi 

It is known [5] that the non-fading LPF for the above example is equal to 2/3. From the graph, we observe 
that for smaller values of p, F-LPF for above hexagon network with fading is greater than LPF with out fading 
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Example A 



Example B 



Fig. 2. Interference graphs for the two example networks 




0.3 0.4 0.5 0.6 0.7 

Probability of a link being 'ON' (p) 



Fig. 3. Bounds on the fading local pooling factor for the Hexagon network 



structure. As p tends to zero, the fraction of time network remains a cycle also tends to be small and it is known 
that GMS is optimal for tree networks. Therefore, it fits well with intuition to see that fading enhances the F-LPF 
for graphs with cycles. 

Example B: A network where fading structure worsens the relative performance of GMS: Consider the graph 
with 3 links a, b, c as shown above. The interference sets for each link is: X a = {6},Z& = {a, c}and I c = {b}. We 
assume each link is either in state 'ON'(l) or 'OFF'(O). So the global channel state '110' denotes that link a and 
b are in 'ON' state and link c is in 'OFF' state. The fading structure is defined as follows: 7r('110') = 7r('011') = 
vrClll') = 1/3. 

For each global channel state, the possible maximal independent sets are as follows: 
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M ab , 



abc 



and 



M, 



bc,abc 



(i o\ 

1 

(° °) 

1 



and 



M, 



abc 



1 

Any vector that belongs to <&{{abc}) can be represented as follows, 

$ = ^M ab [a !-«]' + ^M 6c [/3 1 - p]' + ^M abc [ 7 1 - 7 ]'. 



(10) 



Let 0i be obtained using (a,/3, 7) = (1,0,0) and 2 be obtained using (a,j3,j) = (1/2,1/2,3/4). Evaluating 
the above expression using the above values, we have <\>\ = |[1 1 1]' and 02 = ^[11 1]'- Observing the fact that 
1 02 = 01. using Theorem we have that local pooling factor for the wireless network with the above fading 
structure is less than or equal to |. But, it is known that the local pooling factor of GMS for tree networks (with 
no fading) is 1. 

This result though sounds counter-intuitive, stems from the fact that we allow the fading to be arbitrary. Thus 
fading can act as adversary and as demonstrated, can degrade the performance of GMS algorithm. 



C. Characterization in terms of Interference degree 

So far, we have characterized the performance of GMS through a single scaling factor of the entire throughput 
region. Note that each fading state J induces a network defined on the set of edges that are in 'ON' state and 
GMS can stabilize the network if arrivals are inside the region c*(J)Aj. It is natural to ask for the fading scenario, 
i.e. network with distribution vr(J), if GMS could stabilize the region ^ j7r(J)a*(J)Aj? We answer the above 
question in two parts. 

In the first part, we show the interesting result that GMS cannot stabilize the above averaged region. In other 
words, there exists an arrival process with rate outside the region Af(x) for x(J) = cr*(J) (standard LPF) that can 
make the network unstable under GMS algorithm. We illustrate this using a simple example described below. 
Counter Example: Consider the network with 3 nodes as in Example B. Note that the standard LPF CO for all the 
three fading states is 1. Thus the region Af(a*(J)) is exactly same as the actual throughput region A/. However, 



12 

we have shown earlier that F-LPF is strictly less than 0.8. Thus there exists an arrival process with rates outside 
the region 0.8A^ that cannot be stabilized by the greedy maximal schedule. 

Given the previous negative result, in the second part we show that GMS can stabilize the region A.f( d )j\ )- 
Note that this region is strictly inside the region Aj(x) with x(J) = cr*(J). More formally, our result is as follows: 

Theorem 4. Under a given network topology and channel state distribution with Assumption Al on the arrivals and 
fading channels, GMS can stabilize the network if the arrival rates are inside the region A.f(x), where x(S) = d , s y 

Implications: The above theorem provides an elegant characterization of the rate region that can be stabilizable 
by the GMS algorithm. Also, we find that that the above region is not a subset of the achievable region stated in 
Theorem lb (i.e o-* G (ir)hf). We illustrate the above observation through a simple example described below. 

Consider the wireless network with 3 nodes and fading distribution similar to example B. Note that the interference 
degree for fading state '110' is d/('110') = 1, for state 'Oil' is dj('Oll') = 1 and for the fading state '111' is 
di(' 111') = 0.5. Any arrival rate vector that belongs to the new region defined using the interference degree can 
be expressed as below, 

A = ^M ab [a l-a}' + ^M bc [(3 1 - (3]' + ^M a6c [ 7 1 - 7 ]', (11) 

where a, /3 and 7 are positive constants that are bounded by 1. Using (a,f3,j) = (0, 1,0), we have that rate vector 
(0, |, 0) is inside the new region characterized by the interference degree. However, note that we have shown the 
F-LPF is upper bounded by | for example B network. Thus, all arrival rates that are inside the region |Aj satisfy 

| and hence rate vector (0, |, 



the constraint that A 2 < I and hence rate vector (0, 1,0) belongs to the new region and not the region characterized 



by F-LPF. 

Proof Discussion: We consider the continuous time model with deterministic arrival and channel state processes. 
We then exhibit a Lyapunov function, sum of squares of queue lengths, whose derivative is strictly less than zero 
under the GMS policy whenever the arrival rate is strictly inside the new region. Therefore, the fluid model is 
stable and thus using the results from |f2] we conclude that the original network model is stable. 

IV. Proofs of Results 

Theorem (1). a) (Upper Bound) Under a given network topology and channel state distribution with Assumption 
Al on the arrivals and fading channels, the efficiency ratio of GMS (7*) is less than or equal to o~q(tt). 

b) (Achiev ability) Under a given network topology and channel state distribution tv with Assumptions Al and 
A2 on the arrivals and fading channels, the efficiency ratio of GMS (j*) is greater than or equal to o~ g (tv). 

Proof: The proof follows the method developed by the authors in /[5]/, fi3§ for the non-fading case; however 
we have extended it to take in to account the fading structure. First, for the converse (to show instability for 
arrivals outside the stability region), we explicitly construct an adversarial channel variations pattern that satisfies 
the time-averages imposed by the fading assumption, and this is used in conjunction with the adversarial arrived 



13 

process. The achievability part is more straightforward — we augment the analysis in /|3]/, /[5]/ to include the fluid 
limit of the channel fading process. We now provide the proof more detail: 
Proof (Theorem 1. a): The result follows from the following general lemma. 

Lemma 1. If there exists a subset of links L(C fC), a positive number a and two vectors {1, V € &(L) such that 
ajl > V, then for arbitrary small e > 0, there exists a traffic pattern with offered load u + e&L and a fading pattern, 
such that system is unstable under greedy maximal schedule. 

Proof (Lemma [l}: The idea of the proof is as follows - we construct a traffic pattern and channel variations 
pattern with offered load v + ee^ and show that under this traffic/channel fading pattern, the queue lengths go to 
infinity under GMS, thus making the system unstable. 

As remarked earlier, this proof technique was introduced in fi5§, where authors only needed to construct adver- 
sarial arrival process that makes the queues in the system to overflow. However, in our setting, we need to account 
for the fading process and construct both arrived and channel fading pattern that makes the network unstable. 

Since v G &(L), there exist vectors w J such that v can be expressed as, 

v=Y,kl(J)(Mj, l w j ). (12) 

JCL 

Fix S > 0, we then find a vector f J in the set of rational numbers, Q, such that \\f J — w J \\ < S. 

Assume packets arrive to a link at beginning of the time slot. Let the queues of all the links in L are empty 

at t = 0. Let Tj be the smallest integer such that for all i, rfTj is an integer. Let tj = rfTj. Also, there exists 

integers ni,n,2, ...n 2 L such that 

I njTj i r\\ s S ni\ 



Es-.SCL n sTs 2 L 

Let us define ttl(J) G Q as follows, 

*l(J) := " JT V (14) 

Lscl n sTs 

Using the rational quantities ttl(J) and r J , we define z7 r as follows, 

^= J2 *L{J)(M J)L r J ). (15) 

J-.JCL 

Consider a total time period of J2.j n jTj- We assume that channel state remains in J state for Tj time slots 
(denoted as a time frame). It is easy to observe that with the above described fading pattern, we achieve the same 
channel state distribution as 7Tl(J) on links of set L. We now describe the arrival pattern for Tj time slots when 
the channel is in state J. 

Assume that all the queue lengths (of links in L) are equal at the beginning of Tj time slots. We now construct 
arrival pattern that keeps the queue lengths of all links in set L equal at the end of Tj time slots under the GMS 
policy. The arrival process is as follows: 
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1) The time frame of Tj slots is further divided in to £/,£<{> ....t\isj\ time slots, where tf = rfTj and \IS J \ 
denotes the number of columns in Mj. 

2) During the tf , i ^ \IS J \ time slots, apply one packet to each link that is 'ON' in the i th column of Mj. For 
the last tjjcji trme slots, apply one packet to each link that is ON in the last column of Mj at the beginning 
of the time slot except for the last one time slot. For the last one time slot, with probability 1 — e we do the 
same as described before and with probability e, we apply two packets to each link that is ON in the last 
column of Mj and 1 packet to rest of links in L. 

Note that the arrival process is modified compared to one proposed in [5] so as to ensure that all queues remain 
equal after Tj time slots. 

It is now easy to see that at the end of Tj time slots, all the queue lengths are equal and increase by 1 with 
probability e. Thus the above arrival and channel variation pattern make the system unstable under GMS schedule. 
We now show that the arrival rate is same as u + ee^. 

Let e*j denote the vector of all zeros except for i th position which is set to one. Let J2j = S jcl f° r tne remaining 
part of the proof. For the constructed adversarial arrival process, the arrival rate is given by the following, 

^adv — iTF7| — ', ( 1(D ) 



Rewriting the above expression in terms of #£,(</), we have that 

Tj 



\IS J \ 

A adv = 5>*(J)( E riMje.) + e( £ ^)* (IV) 

j *=i 

Thus we have, 



Aadv = J2^W M J,Lf J )+e(J2 ^f^K (18) 

J J J 

We choose small enough S so that the arrival rate is strictly less than v + ee^. 

Proof (Theorem 1. b): This proof is a simple extension of that in /f5]/, /fj]/, however modified to include the fluid 
limit arising due to the channel fading process. Thus, we have provided a detailed sketch and refer to /O/, /El/ for 
full details. 

We consider the fluid limit of the queuing process and we provide a Lyapunov function and show negative drift 
under GMS schedule whenever arrival rate A € ((Tq(tv) — e)Af. 

Consider a sequence of systems -Q n (nt) (scaled in time and space by a factor of n), where Q n (.) denotes the 
queue lengths of original system, satisfying ^ Qf(0) < n at time t = 0. Let us index the sequence of systems 
by n = {1,2, ....}. We apply the same arrival processes to all the above defined systems (i.e A n {.) = A(.)) and 
assume that queues are served according to greedy maximal schedule. Let A n (t) and D n (t) denote the cumulative 
arrival and departure process of system n up to time t. 

Using the results from 0, it can be shown that the sequence of processes (Q n (.) , A n (.) , D n (.)) as n — > oo 
converges to a fluid limit almost surely along a subsequence {?ifc} hi the topology of uniform convergence over 
compact sets, 
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— A?*(n k t) -> A,t, (19) 

rik 

—D^(n k t) -► J>(J)( / ri(s)ds), (20) 

n fc j Jo 

—QT(n k t) -> «(*)• (21) 



Also, the fluid limits (qi(t),/j,f(t)) satisfy the following equality: 



<fi(<) = <ft(0) + *J*-y>(-7)( / /jfto&O- (22) 

Moreover, fluid limits are absolutely continuous, and at regular times t (i.e., those points in time where the 
derivatives exist) we have the following condition satisfied: 

d , . f A, - mi(*) if «(*) > 

dt ( (A, -#(*))+ if 9/(0 = 0, 

where ^(i) = Sj 7r (^)/ i f (*) satisfies the GMS properties. Let Lo denote the set of links with the longest queues 
at time t, 

L (t) = {i € K\qi(t) = max jeK qj(t)} (23) 

Let L(t) denote the set of links with the largest derivative of queue length among the links in Lo(t), 

L{t) = {i£ L (t)\- qi (t) = max i€Lo{t) - qi (t)} (24) 

Lemma 2. Under the greedy maximal schedule, the service rate satisfies /Z(t)|x,(t) £ ®{L(t)), where u\l denotes 
the projection of vector on u on to set of links L. 

The proof of the above lemma is similar to one in J3], and is presented in appendix. The idea, roughly is 
that, queues in the set L(t) will remain the longest for small enough amount of time past t and GMS picks the 
maximal schedule restricted to links in L(t) that are in 'ON' state. 

Since the arrival rates are strictly with in a* L {ir)Kf, there exists a service vector V € &(L) and v < o-* L {it)hf 
such that X(L) < v, where \(L) is projection of arrival vector on to the set L. Given any two vectors in set &(L), 
note that one vector never dominates the other one in all the dimensions by a factor more than a* L {ir). Therefore 
we have that ^rnaXjg^mg^t) is strictly negative when ever maxqj(t) > 0. 

Let V(t) = max qi(t) denote the Lyapunov function used for the fluid system. Since we have a negative drift for 
the Lyapunov function, using the results from J2J, we have that fluid system is stable (i.e there exists to > such 
that qi(t) = OVi > to)- Therefore from |f2], we have that the queues in the original queuing system are stable. ■ 
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Theorem (2). For every J C K, and any (/2j, vj, Hj) such that jlj, uj € CH(Mj), uj < Hjflj, we have that 

a G (7v) < maxi -^== , 

where [Xj(l) = if I £ J. 

Proof- 
Since (iTj,vj,Hj) satisfy the inequality, 

vj < Hjflj (25) 

Summing over all subsets with positive scaling constants it (J), 

Y, *(J)vj(l) < E 7r ( J ) ( H jM0) (26) 

J J 

Using the maximum constant over all the inequalities, we have the following, 



By observing the fact that (Ylj^iJ^j, Sj 7r ( t 0$j) De l° n g to the "£(£)> we have the result. ■ 

Theorem (3). 

., w SjcL^(J)n(Mj) 

where n(M) = minj ^- Mij,N(M) = rnaxj ^- My am/ 7T£, denotes the marginal distribution on set of links L 
induced by it. 

Proof: We first state a lemma that describes the dual problem that finds the fading Local Pooling Factor as 
the optimal solution. The dual characterization of Local Pooling Factor was presented previously in ||3], ID. We 
now provide such characterization for F-LPF in Lemma [3] by generalizing the arguments in J9J. In particular, the 
multiple global channel states due to fading each induce a different constraint - combining all of these appropriately 
while satisfying the long-term average fractions {ttl(J)} results in a maxmin problem, as detailed below. This 
result is used to derive the lower bound. 

Lemma 3. The following optimization problem characterizes o~* l {tv) : 

o~l(iv) = max 2_, ^l(J)o,(J) 
j-.jcl 

s.t : x'M JyL > a(J)e' VJCL 
x'M JtL < b(J)e VJCL 
5>x(J)6(J) = 1 

JCL 
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Proof: Consider the definition of ct* l (it) in I®. The corresponding optimization problem is given by: 

inf a 

s.t : a Y, 7T L (J)Mj tL a(J) > ^ vr L (J)M JiL /3(J) 

JCL JCL 

||a(J)|| = l V JCL 

||/3(J)|| = 1 V JCL 
a{J)J(J) >0 

where |.| is defined as the sum of all the elements of the vector. Let us define a new variable j(J) = aa(J). 
Thus, we have: 

inf a 

s.t : Y, vr L (J)M JiL (/3(J) - f(J)) < 

JCL 

||7(J)||=ct V JCL 

||/3(J)|| = 1 V JCL 
l(J)J(J)>0 

For the above LP, let (x, {y(J)}, {z(J)}) denote the dual variables associated with the constraints. The dual is 
given by 

max min er+ 

£,{v(J)},{z(J)} <r,5(J),/3(J) 

L \ISj\ 

J>(j> L (J)[£M^(/3/- 7 /)]) + 

i=l JCL j=l 



Yy(J)(l(jye-a) + 

JCL 

5>(J)(/3(J)'e-l) 



JCL 

s.t:7(J),/3(J)>0 
Rewriting the above dual optimization problem, we have 



,. ,™, m min . - XX J ) + ^ - X^( J ))+ 

,{v(J)},{z(J)} ct,3{J),P(J) j j 

\ISj\ L 

^/3/[vr L (J)^x J M^ + z(J)] + 

j=l i=l 

\ISj\ L 

j2-ft[-KL(J)J2 x i M £ + yW} 

j=l i=l 

s.t:7(J),/3(J)>0 
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Equivalently, the above program can be reduced to 

E -*( J ) 



max 

j-.jcl 



s.t : ir L {J)x'M J)L + z(J)e' > VJCL 
- tvl(J)x'Mj, l + y(J)e > VJ C L 

E v( J ) = l 

JCL 

Denoting ~3j\ by a(J) and ^4-4 by b( J) we have the desired result. 

■ 
From the above Lemma [3] we have that a"2(7r) is equal to, 

max > iri(J)a(J) 

s.t : i'Mj,i > 0(7)6* VJ C L 
x'M/ jL < 6(J)e' VJCL 
J> L (J)6(J) = 1 

JCL 

Observe that ( ^ TJWTm ) e > T n (j)n(m ) ' 1) * s a van d point in the search space. Substituting the point in the 
above function, we have the desired inequality. ■ 

Corollary 1: a* G {ir) > j^ 

Proof: Observing the fact that n(Mj) > d l G , N(Mj) and using the above lemma, we have the desired 
inequality. ■ 

Theorem (4). Under a given network topology and channel state distribution with Assumption Al on the arrivals and 
fading channels, GMS can stabilize the network if the arrival rates are inside the region Af(x), where x(S) = d , s y 

Proof: We consider a continuous model similar to the one described in the proof of Theorem lb. In this model, 
the queuing system evolves according to the following equation, 

d u , \ h~ Mt) if «(*) > ° 

dt \ (\l-^l(t)) + if qi(t)=0, 

where p,i(t) = Sj 7r (7)A i f (*) satisfies the GMS properties. In the original system with fading channels note that 
the weight of GMS schedule is always greater than d 1„, of the weight of the max-weight schedule where S is 
the set of links that are in 'ON' state. Therefore in the fluid model, we can show that ///(£) satisfies the following 
condition 

2 dl{J) yj&CH{Mj, K )'-r' 



19 



Let us the consider the following Lyapunov function, 

I 
Taking the derivate of the Lyapunov function, we have that 

WHS^ifW-lnW). (30) 

l 

Using the GMS properties of m (t) , we have 

TA(^))<(2^«(*)Ai-X)jTn 7r ( J )- £3* X«^ J WJ (3D 

As A is assumed to lie inside the region Aj(x), there exists ffj € CH(Mj i )c) such that 



^<E^j)-( J )^(0- 



(32) 



Using the above inequality, we have that 



V(q(t)) <(2$>(*)E TTn^^W "Ett^ j ) - £S* X*^ J W) (33) 

Thus from the above inequality we have that V(q(t)) < whenever q(t) > 0. 

We can now use the results from ||2] to argue that the original system is stable under the assumed arrival process 
as the fluid model is stable. ■ 

V. Extensions to Multiple Fading States 

We now extend our results for 'ON/OFF' channels to channel models where each link capacity is time- varying 

and takes values from a finite state space. Let us denote the set of values in the state space by {0, c±, C2, , c m }. 

The global state GS(t) of the system now refers to the exact channel state of each link. Let n(X±,X2, ...,Xk) 
denote the fraction of time the network is in global channel state (X\,X2,Xs, ....Xk)- Let us denote the state 
(X 1 ,X 2 ,X 3 ,....,X K )byX. 

Let Mx denote the matrix consisting of K rows one for each link. Each column now represents a possible 
maximal independent set on the set of links with non-zero channel states. For a given column, the entries of a 
given row is set to zero if link I (corresponding to row) does not belong to independent set, or is set to equal to 
channel value X\ if it belongs to independent set. For example, consider the Interference graph in Figure J] with 
each link taking 3 channel states {0, 1, 2}. Then Mn ,2,1.0) i s given by, 



JW(1,2,1,0) 



(1 


o\ 





2 


1 





v° 


°; 
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The throughput region Ay for the above general network model with fading pattern 7r(X) is given by: 

A 9 f = {X:X>0 , A<^ vr(X)r? x where 

x 
r/ x € CH(Mx)}. 

We now define the F-LPF for a set of links L as follows: 

cr* (tt) = inf{<j : 3 <&, <f 2 G $ 9 (£) such that a& > 2 }, (34) 

where, 

*ff(L) = {0 : <£= ^ 7r(X)fjk where 7? x € C«(M Xl )}, (35) 

x 

Xl is constructed from X by setting the values of links that do not belong to set L in X to zero. 

Theorem 1 can be shown to hold for the general model with the above modified definition of F-LPF. The proof 
of Theorem 1 for the 'ON/OFF' channels can be easily modified to above system with general channels and is 
therefore omitted. 

VI. Conclusion & Discussion 

In this paper, we studied the problem of scheduling in wireless networks with interference constraints where the 
capacity of links changes over time. We have analyzed the performance of a well-known algorithm, Greedy-Maximal 
Scheduling (GMS), to the case of general wireless networks with fading structure. We defined Fading-Local pooling 
factor for graphs with fading and showed that it characterizes the fraction of throughput that can be achieved by 
GMS. We have derived useful yet easily computable bounds on F-LPF through alternate formulations. 

By analyzing F-LPF, we have studied the effect of fading on the performance of GMS. It is a priori not clear 
whether fading can enhance/degrade the relative performance of GMS. In this work, we have showed that fading 
can in fact exhibit both behaviors through two simple examples, one in which fading increases the efficiency ratio 
of GMS and other in which fading decreases the efficiency ratio as compared to non-fading case. 
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VII. Appendix 

Proof: (Lemma |2]) 
The proof is similar to the one presented in ||5] however taking in to account the channel fading. From the 
definition of set Lq{€) in Eqn (|23T ). there exists ei > such that 

Qi(t) >Qj(t)+ei Vi £L (t) and j eK\L (t). 
Using the continuous property of qi(t), we further have that, there exists e 2 > 0, 8\ > such that 



min qAt + 5) > max qAt + 8) + eoV£ £ [0, £il. 
i€Lo(t) jeK\L (t) 

Since L(t) is contained inside Lo(t), we have that, there exists e 2 > 0, 8\ > such that 



min qAt + 5) > max qAt + 5) + e 2 V5 G [0, 81]. (36) 

ieL(t) jeK\L (t) 



Also, from the definition of set Lit) in Eqn d24b . there exists €3 > such that 



j t qi(t) > j t Qj{t) + £3 ViG L(t) and j G L (t) \ L{t). 
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Further, using the definition of derivative ^q(t) « — + ' s ~ 9 ^ - , there exists e^> 0,62 > such that the following 
holds. For all i G L(i) and j € Lo(£) \ L(t), we have 

^ + ^-^ ) > ^ + ^-^ ) +g4 v^(o^ 2 ] 

d d 

Using the fact that queues qi(t) in set L (t) are equal, the above inequality can be rewritten as follows. For all 

% G L(t) and j G Lo(t) \ L(t), we have 

qi(t + 5) q 3 (t + 5) yXtz , nx , 
7 > -^— ? he 4 Ve) € (0, d 2 J- 



Thus we have, 



min qi(t + 6) > max a At + 6) + e 5 V£ G (O.fol. (37) 

»6L(t) jeL (t)\L(t) JK 



From the inequalities (l36T l and (l37l . we have the following inequality, there exists 60,63 > such that for all 

6 G [^0,53] we have 

min qAn(t + 6)) > max qAnit + <5)) + es. (38) 

iei(t) jelc\L(t) 

From the definition of fluid limit qi(t), there exists uq large enough such that Vra > no and<5 G [6q, 63], we have 

that 

min Qi(n(t + 6)) > max QAn(t + 6)) + ne 7 . (39) 

i€i(t) jelC\L(i) 

The above inequality ensures that the links in the set L(t) have larger queue lengths compared to other links in 

the network for all the time slots in [n(t + 60), n(t + 63)]. Therefore, depending up on global channel state GS(t), 

at each time slot r G [n(t + 6o),n(t + 63)], GMS schedule picks a maximal schedule from the set of links L(t) 

that are in 'ON' state. Let Zf(r) denote the scheduling decision picked by the GMS algorithm for link / at time 

slot r. We thus have 

Z n (T)\ L[T) eM GS[T)nmMt) . (40) 

Computing the total service provided by the GMS algorithm in time slots \n{t + 60), n(t + 63)], we have 



rnt+nS 3 

Df{nt + n6 3 ) - Df{nt + n6 ) = / Z?{t)<1t. 

Jnt+nSn 



Let us denote the quantity ' (n — n 3 A — pp — ^-^ by ftf(t). From the above equality, we have that /£"(£) |_L(t) G 
<fr(X(£)). As 6q can be made arbitrarily small, we have the result. 



